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ABSTRACT 

We present a modified scenario of gamma-ray emission from pulsars within the 
framework of polar cap models. Our model incorporates possible acceleration of 
electron-positron pairs created in magnetospheres, and their subsequent contribution 
to gamma-ray luminosity L 1 . It also reproduces the empirical trend in L 7 for seven 
pulsars detected with Compton Gamma Ray Observatory (CGRO ) experiments. At 
the same time it avoids basic difficulties (Nel et al. 1996, Arons 1996) faced by theo- 
retical models when confronted with observational constraints. 

We show that the classical and millisecond pulsars form two distinct branches 
in the L 7 - L s( j diagram (where L s( j is the spin-down luminosity). In particular, we 
explain why the millisecond pulsar J0437-4715 has not been detected with any of the 
CGRO instruments despite its very high position in the ranking list of spin-down 
fluxes (i.e. L s d/D 2 , where D is a distance). The gamma-ray luminosity predicted for 
this particular object is about one order of magnitude below the upper limit set by 
EGRET. 

Key words: gamma-rays: theory, observations - pulsars: general 



1 INTRODUCTION 

Numerous models of pulsars (see Michel 1991 for the re- 
view) proposed over the last three decades tried to make 
specific predictions about emission of gamma-rays and X- 
rays. Gamma-rays are particularly important as a direct 
signature of basic non-thermal processes in pulsar magneto- 
spheres, and potentially should help to discriminate among 
different models. Interpreting gamma-rays should also be 
less ambiguous compared to X-rays. In the latter case, es- 
pecially for objects younger than 10 6 yr, contributions from 
initial cooling, internal friction, etc. of unknown magnitude 
may dominate the total X-ray emission. 

There are seven positive detections of pulsars by 
CGRO (see Table 1), i.e. less than one per cent of all pulsars 
known to date. In all cases the sources had been identified 
by virtue of gamma-ray flux modulations with previously 
known P. Crab and Vela are the only pulsars seen by three 
of CGRO detectors. 

The EGRET data became so far the only ground for 
testing theoretical models. The latest critical review of three 
models of gamma-ray emission (polar cap models by Hard- 
ing 1981, Dermer & Sturner 1994, and outer gap model of 
Yadigaroglu & Romani 1995) comes from Nel et al.(1996). 
They conclude that none of the models fits observations sat- 
isfactorily for two major reasons. First, the confrontation 



of polar cap models (Harding 1981, and Dermer 
1994) with the observations leads to a relation 



L 7 (observ.) tx L"(model), 



Sturner 



(1) 



with a ~ 0.6 and 0.5, respectively, instead of a = 1. Sec- 
ond, there have always been some cases amongst ~ 350 
EGRET upper limits, apparently contradicting predictions 
of L 7 made by the models. The troublesome limits come 
usually from the pulsars B1509-58, B1046-58, B0656+14, 
B1929+10, B0950+08, as well as from the millisecond ob- 
ject J0437-4715. 

Of these two problems raised by Nel et al.(1996), the 
latter is more severe in our opinion. The former problem 
may be solved to some degree by updating parameters used 
for two objects at the low- luminosity domain, i.e. B1055- 
52 and Geminga. After Thompson et al.(1994), Nel et al. 
used for B1055-52 the spectral index 7 = 1.18 determined 
by Fierro et al. (1993) from the first three viewing periods of 
EGRET . However, a substantially higher value, 7 = 1.59, 
based on the data from 10 viewing periods became recently 
available (Fierro 1995) . This value of the spectral slope re- 
duces L-y of B1055-52 by a factor of ~ 4. (We shall discuss 
other consequences of the steeper spectral slope for B 1055-52 
in Section 3) . Further reduction of L 7 in the case of B 1055-52 
is possible by lowering distance D to the source. The argu- 
ment for lowering the distance (the usually assumed value 



© 0000 RAS 



2 B. Rudak & J. Dyks 

Table 1. Gamma-ray luminosities for pulsars detected with CGRO (in log of [ erg s — 1 ]). Beaming angle of emission Q 7 = 1 sr was 
assumed. 



PSR 




log L sd 


EGRET 


COMPTEL 


OSSE 


Refs. 


'totaF 


B0531+21 


Crab 


38.65 


34.6 


35.0 


35.1 


1,2,3 


35.42 


B1509-58 




37.25 






34.4 


-,-,3 


34.40 


B0833-45 


Vela 


36.84 


34.2 


33.5 


31.3 


1,2,3 


34.28 


B1951+32 




36.57 


34.1 


34.1 




1,4,- 


34.40 


B1706-44 




36.53 


34.4 






1,-r 


34.40 


J0633+1746 


Geminga* 


34.51 


32.9 


+ 




1,5,- 


32.90 


B1055-52 




34.48 


33.4 






6,-r 


33.40 



*at the distance of 157 pc. 
Refs. to flux values: 

1) Ramanamurthy et al.(1995), 2) Carramifiana et al.(1995), 3) Schroeder et al.(1995), 
4) Kuiper et al.(1996a), 5) Kuiper et al.(1996b), 6) Fierro (1995) 



is 1.53 kpc, after the model of Taylor & Cordes 1993) may 
come from ROSAT PSPC observations. If the bulk of the 
(presumably) thermal X-ray emission from B1055-52 is due 
to initial cooling then the inferred radius for the neutron 
star at D = 1.53 kpc exceeds 30 km (Ogelman 1995) - a 
value hardly acceptable from the theoretical point of view. 
In the case of Geminga, its inferred luminosity drops by a 
factor of ~ 2.8 relative to the value used by Nel et al. due to 
a newly determined HST parallax distance of 157 pc (Car- 
aveo et al. 1996), instead of D = 250 pc. Fig. [j] presents 
the inferred EGRET luminosities versus model predictions 
of Harding (1981), with Geminga, and B1055-52 points cor- 
rected for D = 157 pc, and 7 = 1.59, respectively. Generally, 
the agreement between theory and observations looks quite 
satisfactory, even though the Crab pulsar visibly deviates 
from the trend predicted by Harding (1981). 

Arons (1996) presented a simple argument against all 
models which offer a functional relation for L 7 roughly 
similar to L 7 oc B/P 2 , by plotting the voltage (~ B/P 2 ) 
available for particles, against the gamma-ray luminosity 
L 7 (> 100 MeV) for six EGRET pulsars. An extrapolation 
of the resulting trend towards low values of B/P 2 leads to 
L 1 exceeding L s d below ~ 10 14 Volts. Since L 1 must never 
exceed L B a, this special value of voltage defines then 'an 
empirical gamma-ray death line' on P — P diagram (or on 
Voltage — P diagram, as originally presented by Arons). It is 
hard to reconcile such a line with the observed death line for 
radio emission, which corresponds to ~ 10 12 Volts. The effect 
was known already earlier, when polar cap models of Buc- 
cheri et al.(1978) and Harding (1981) had been introduced. 
In the latter model, for a pulsar with P ~ 10~ 15 s s~\ its 
L 7 would reach L s d at a characteristic age equal to 3 x 10 7 
years (Harding 1981). 

Below we propose a polar cap model which is free of 
the two problems discussed above. This model reproduces 
gamma-ray luminosities inferred for seven observed pulsars. 
At the same time it avoids the problem of the empirical 
gamma-ray death line of Arons (1996), and it relaxes the 
upper limits' constraints of Nel et al. (1996), especially for 
old classical pulsars, and millisecond pulsars. In Section 2 
we start with recalling the model for total power contained 
in outflowing particles, which refers directly to the relation 
analysed by Arons (1996). Then we present arguments for 
summing all available CGRO data in order to get a better 




32 34 36 

22. + log (L s / 3 7t° 15 ) [in erg/s] 

Figure 1. Gamma-ray luminosity above 100 MeV for 6 
EGRET sources is plotted against predictions of the polar cap 
model of Harding (1981). The normalization factor is arbitrary; 
t = P/2P is a characteristic age expressed in [yr]. (Note: The 
Crab pulsar is the first dot from the right.) Bars indicate uncer- 
tainties in L 7 arising from uncertainties in distance and flux. In 
the case of Vela (next to Crab) the lower distance limit was ex- 
tended down to 200 pc to conform to recent estimates based on X- 
ray and gamma-ray observations of the Vela SNR (400±200pc af- 
ter Aschenbach et al.1995 and Si 350 pc after Oberlack et al.1994, 
respectively) . 

start for a modified polar cap model rather than using the 
EGRET data alone. Section 3 contains a description of our 
model, and its reference to existing information on gamma- 
rays from pulsars. Summary and comments are in Section 
4; it contains also the ranking of pulsars with the highest 
gamma-ray fluxes resulting from our model. 
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2 SIMPLE MODELS VS. CGRO DATA 

Though EGRET has been far more successful than other 
CGRO instruments in detecting pulsars, this does not mean 
that pulsars' gamma-ray emission above 100 MeV domi- 
nates energetically over gamma-ray emission from lower en- 
ergy bands. In the case of the Crab pulsar most of the en- 
ergy output occurs within the COMPTEL and OSSE energy 
ranges (e.g. Fierro 1995). In the extreme case of B1509-58 
there is only OSSE detection. Early positive reports from 
the COMPTEL team (Carraminana et al. 1995) have not 
been confirmed (Kuiper 1996) , and EGRET put only upper 
limits for the source (Thompson et al. 1994). Luminosities 
listed in Table 1 were inferred from phase-averaged fluxes as- 
suming beaming solid angle of gamma-ray emission f2 7 equal 
1 steradian. For 3 objects (Crab, Vela, and B1951+32) de- 
tected with more than one instrument, we can also construct 
a sum of inferred luminosities (ignoring possible changes of 
a beaming angle with energy range), to get an estimate of 
'total' gamma-ray luminosity. For the remaining 4 pulsars 
we will use their luminosities in the EGRET energy range 
as 'total'. The values of 'total' gamma-ray luminosities are 
shown in the last colum of Table 1. 

Let us now compare these 'total' gamma-ray luminosi- 
ties with the simplest possible phenomenological polar cap 
model. According to this model, gamma-ray luminosity L 7 is 
proportional to a power contained in outflowing primary 
electrons I/ part icics, which in turn is proportional to a prod- 
uct of primary electron energy Eo, a surface area of canon- 
ical polar cap A pc ~ 7r.Rp c oc l/P, and a Goldreich- Julian 
flux nGJ oc B/P (Goldreich & Julian 1969) of outflowing 
primary electrons: 

L-, OC Lparticlcs OC Eo ■ flGJ ■ A pc . (2) 

Assuming Eo = const, for all objects one obtains 

Lry OC Lp ar ti c l cs OC B/P 2 oc L\' d 2 . (3) 

(Note: The model for L part icics from equation (^) was ac- 
tually a starting point in the work of Harding 1981, who 
assumed Eo = 10 13 eV and whose objective was to find a 
prescription for L 7 (> 100 MeV). We shall return to this 
point in the next section.) 

Figure ^ shows how the relation of equation (^) com- 
pares with observations. The normalization factor (C = 
10 16 ) has been chosen to obtain the best 'by eye' fit. The 
overall agreement looks quite impressive. There is a sub- 
stantial improvement for Crab comparing to Figure |l| due to 
significant contributions from OSSE and COMPTEL data. 
Moreover, B 1509-58 adds up smoothly to six EGRET pul- 
sars. The same functional relation, but for EGRET points 
only, has been considered by Arons (1996) (see previous sec- 
tion). 

The relation L 1 = 10 16 Lgj 2 [ erg s" 1 ] presented in 
Fig.^, cannot hold for all radio pulsars. It leads formally to 
L-y reaching L a< j at 10 32 erg s _1 (which corresponds to the 
Arons' empirical gamma-ray death line), whereas pulsars are 
observed down to L s( j — 10 30 erg s . Clearly, pulsar mod- 
els which predict L 7 as a simple combination of B and P, 
require some revision. 




32 34 

16. + log L sd ' [in erg/s] 



Figure 2. 'Total' gamma-ray luminosity L 7 for all seven 
CGRO sources, inferred from EGRET , COMPTEL , and 
OSSE observations (wherever available - see Table 1.) are com- 
pared with a simple model (equation g|) discussed in Section 2. 
(Note: B1951+32 and B1706-44 practically coincide. The dot ly- 
ing at the upper right part of the diagonal represents the Crab 
pulsar.) 

3 HOW DO ELECTRON-POSITRON PAIRS 
CONTRIBUTE TO GAMMA-RAYS 

According to the model of Daugherty & Harding 1982 
(DH82) primary electrons are accelerated along open mag- 
netic field lines to high energies (~ 10 13 eV) due to rotation- 
induced electric field. The model assumes a dipolar struc- 
ture of the magnetic field. Curvature photons emitted by 
primary electrons are absorbed by magnetic field with sub- 
sequent creation of electron-positron pairs (Sturrock pairs). 
These pairs cool off instantly via synchrotron radiation. Syn- 
chrotron photons may lead to further pair creation. Elec- 
tromagnetic cascades propagating in pulsar's magnetic field 
may be very rich, with several subsequent generations of 
pairs and photons. 

Numerical treatment of electromagnetic cascades initi- 
ated by primary electrons above polar cap, and propagating 
across the magnetosphere (DH82) does not include effects of 
possible acceleration of Sturrock pairs. The only contribu- 
tion from pairs to gamma-rays taken into account is due to 
synchrotron emission of created pairs. The pairs themselves 
do not accelerate, and subsequently - do not contribute to 
the curvature radiation. Such simplification is usually jus- 
tified by arguing that an appearence of conductive plasma 
above some height effectively leads to a screening of elec- 
tric field parallel to local magnetic field lines. If, however, 
the density of created pairs is lower than the local corota- 
tion plasma density, the electrons from pairs will be subject 
to further acceleration (whereas positrons will be deceler- 
ated; eventually some of them will be stopped and reversed 
towards the stellar surface). In the context of polar cap mod- 
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els developed by Daugherty & Harding (1982, 1994, 1996) 
it became clear that potential contribution to gamma-rays 
from pairs might be necessary to account for the observed 
gamma-ray fluxes . 

If these secondary particles are indeed subject to effec- 
tive acceleration at significant altitudes above polar cap sur- 
face, e.g. at heights of several NS radii (Daugherty & Hard- 
ing 1996), the resulting beaming angles of gamma-ray emis- 
sion fi 7 will be wider than those measured at the polar cap 
surface. The requirement for "nearly aligned rotators" (Der- 
mer & Sturner 1994, Daugherty & Harding 1994), necessary 
to explain large duty cycles of gamma-ray emission, might 
be then relaxed. [Note: Relaxing the assumption about small 
inclination angles between the rotation and magnetic axes 
may be inevitable on observational grounds. Inclination esti- 
mates carried out by Lyne & Manchester (1988), and Rankin 
(1990) show that in many pulsars inclination angles are large 
indeed. However, mutual comparison of these results shows 
that they are in agreement only for small inclination angles 
(Js 40deg), and there is no correlation if either estimate is 
larger than this value (Miller & Hamilton 1993).] Hereafter 
we will assume that for all pulsars f2 7 = lsr (corresponding 
to opening angles of ~ 30 degrees). 

Suppose, that the secondary particles, which we assume 
to be created as described in the model of DH82, do partic- 
ipate in the gamma-ray production similarly as the primary 
electrons do. In the spirit of equations Q and (j3|), we pro- 
pose a following prescription for L 7 : 

L 7 = C-n± -E± -Lli 2 , (4) 

where n± is a number of created pairs e ± per primary elec- 
tron, and E± is a characteristic energy attained by particles 
due to acceleration. The normalization constant C will be 
determined by fitting the observations. We will assume that 
E± achieved by secondary particles is similar to the energy 
attained by primary electrons Eo, i.e. E± ~ Eo. 

We begin with making a choice for the value of Eo, 
since it is this parameter which, along with B and P, will 
determine the number of pairs n± created per primary elec- 
tron. The analytical fit of Harding (1981) for E 7 above 100 
MeV, which gained so much popularity in testing her po- 
lar cap model against EGRET observations was obtained 
from numerical simulations performed for a fixed value of 
primary electron energy Eo = 10 13 eV. Originally, this en- 
ergy has been chosen to make the best spectral fits above 
100 MeV for COS-B data of Crab and Vela. However, such 
assumption about Eo is not possible throughout the entire 
pulsar's lifetime because the energy of outgoing electrons is 
subject to twofold limitation (e.g. Sturrock 1971): there is 
an absolute upper limit 

E w = 1.2 x 10 7 B 12 P~ 2 [ MeV] (5) 

(B12 = B/10 12 G, period P is in [s]) set by potential drop 
across the polar cap, and a maximum value 

Smax = 4.6 x 10 7 BI(, 4 P- 1/8 [ MeV] (6) 

set up by curvature cooling in purely dipolar magnetic field. 
The energy Eo of primary electrons must conform to one 
of these two limits (whichever comes first) as the pulsar is 
slowing down, since both Ew and E max decrease as P in- 
creases. In consequence, the prescription for L 7 (> 100 MeV) 



of Harding (1981) should not be treated as accurate every- 
where. 

Moreover, accelerating electrons should cross threshold 
energy 

E min = 1.2 x 10 7 B~ 2 1/3 P 1/3 [ MeV] (7) 

required for pair creation. Whenever Ew falls below E m i n , 
a classical pulsar crosses the well known deathline, and en- 
ters 'a graveyard for pulsars'. Another deathline occurs for 
millisecond pulsars even earlier, when E max falls below E m i n 
(Rudak & Ritter 1994). 

An ultrarelativistic primary electron, sliding along 
curved magnetic field line, emits curvature photons which 
in turn may be converted into e pairs. The total number 
of pairs created per electron depends on the local component 
of magnetic field B perpendicular to the direction of propa- 
gation of the photon, on the energy of the electron Eo, and 
on the curvature of magnetic field lines. As the primary elec- 
tron accelerates, its energy Eo will cross the threshold value 
E m i n , thus triggering creation of a first pair. The number of 
created pairs quickly increases as Eo exceeds -Emin- At some 
point it becomes high enough to make further acceleration of 
electrons less effective due to screening effects. Instead of fix- 
ing Eo at some specific value we assume thus that electrons 
may be accelerated up to an energy Eo satysfying following 
condition throughout pulsar's entire lifetime: 

Eo = min{£ ■ E mm , -Ew, E max }, (8) 

where ( > 1. 

The best choice for the value of the parameter £ was 
made a posteriori, to reproduce the empirical trend of L 7 for 
the seven CGRO pulsars with similar accuracy as equation 
(|) does in Fig. 2. We found that the range 2 £ Ss 5 fulfills 
this requirement. All results presented below are for £ = 2.5. 
For -E m in we preferred to take the numerically obtained val- 
ues whenever they differred from the analytical approxima- 
tion (the analytical formulae for Ew, -Emax, and -E m i n are 
taken from Rudak & Ritter 1994). We found that the for- 
mer are consistently smaller by a factor of ~ 1.5 (in most 
cases) than the latter. 

The developement of cascades was followed by means 
of numerical simulations described in DH82. Calculations of 
number of pairs n± were performed with numerical simula- 
tions after choosing P and B, and setting the primary elec- 
tron energy Eo according to equation (^). The normalization 
constant in our model of L 7 (equation ^J) was determined 
by fitting numerical results to the seven detections (the last 
column of Table 1.). Then we calculated two evolutionary 
tracks in the L 7 - L s( j space for representatives of the 
classical pulsars (with typical magnetic field strength B ~ 
10 12 G), as well as of the millisecond pulsars (B ~ 10 9 G). 
Both tracks are shown in Fig.[] as solid curves. The upper 
curve (B ~ 10 12 G) starts at L s( j — 10 39 erg s _1 , nearby 
Crab, and down to L s( j ~ 10 34 erg s _1 it roughly follows 
the dashed line, which depicts the relation L 7 = 10 16 • L^ 2 
from Section 2. As L s( j decreases, our exemplary classical 
pulsar enters a region where proximity to pulsar's death line 
becomes important. The number of created pairs n± declines 
constantly as the pulsar slows down, and it starts to decrease 
dramatically when Ew falls below £E m i n , affecting thus Eo 
in equation (js|). 

At the point where -Ew = E m i n , the creation of pairs 
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Figure 3. 'Total' gamma-ray luminosity is plotted against spin- 
down luminosity L a< j for seven CGRO pulsars (filled dots). Open 
triangles are the EGRET upper limits as given by Nel et al.1996 
for 350 objects, including seven millisecond pulsars. Filled trian- 
gle indicates position of J0437-4715 (after Fierro et al.1995). The 
dashed line corresponds to L 7 = 10 16 ■ L^ 2 (see Section 2). The 
upper solid curve shows the evolutionary track of a pulsar with 
B pc = 10 12 G, calculated according to eqs. (^) and (|i|). The lower 
solid curve shows the evolutionary track for a millisecond pul- 
sar with B pc = 10 9 G. The evolutionary tracks end up when the 
objects reach their 'death points' in the L s( j space (marked with 
dotted lines). 



ceases (n± = 0) - the pulsar reaches its death line. The ex- 
emplary millisecond pulsar (B pc = 10 9 G) follows the lower 
solid curve in Fig.^. It starts at L s d ~ 10 37 erg s _1 , cor- 
responding to initial period of one millisecond. Unlike the 
former classical pulsar, it encounters different death line at 
L s d ~ i0 33 ' 5 er g g-i^ due to an equality -E ma x = -Emm- 
The efficiency of pair creation for the millisecond pulsar 
throughout its lifetime, is significantly lower than for the 
classical pulsar. Low strength of magnetic field plays a de- 
cisive role here, and it cannot be compensated for by faster 
rotation (and therefore - by smaller curvature radii avail- 
able). The filled triangle in Fig.[] indicates the upper limit 
for J0437-4715 set from EGRET observations (after Fierro 
et al.1995). In addition, 350 EGRET upper limits (Nel et 
al. 1996), including seven limits for millisecond pulsars, are 
shown for comparison. Apart from J0437-4715, four objects 
with EGRET upper limits only are placed clearly below 
the upper evolutionary track. These are (from right to left): 
B1046-58, B0656+14, B1929+10, and B0950+08. If the pro- 
posed model is correct, these objects should be the best 
candidates for detection in gamma-rays (but not necessar- 
ily in the EGRET energy range). B1951+32 and B1509-58, 
present in the data of Nel et al. (1996), have been replaced 
with detections (filled dots) by EGRET and OSSE , respec- 
tively. 




30 32 34 36 
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Figure 4. 'Total' gamma-ray luminosity is plotted against pre- 
dicted luminosity, calculated according to eqs. fl4) and (pi, for 
seven CGRO pulsars (filled dots). Open triangles denote com- 
bined OSSE , COMPTEL and EGRET upper limits (wherever 
available). The continuous diagonal line corresponds to a perfect 
agreement between predictions and observations. 

The comparison of how our model reproduces L 7 for 
the seven CGRO pulsars, along with combined upper limits 
from EGRET , COMPTEL , and OSSE (from Thompson 
et al.1994, Fierro et al.1995, Schroeder et al.1995, and Car- 
raminana et al.1995), wherever available, is shown in Fig.^J. 
In the case of Geminga and B1509-58, the model overesti- 
mates the observed L 1 by a factor of ~ 2. However, existing 
upper limits from COMPTEL and EGRET , respectively, 
improve the agreement. The upper limit for J0437-4715, 
based on EGRET only, is one order of magnitude above 
the predicted value of L 7 . For four other objects (B 1046-58, 
B0656+14, B1929+10, and B0950+08) the stringent upper 
limits from EGRET differ from the model predictions by 
no more than a factor of ~ 3. Moreover, upper limits from 
COMPTEL on B1046-58 and B1929+10 place them on a safe 
side of the diagonal line of prefect agreement in Fig.^. There 
is no information available about B0950+08 and B0656+14 
from any COMPTEL observations. 



4 SUMMARY 

We have proposed a semi-phenomenological model of 
gamma-ray emission from pulsars, which is based on po- 
lar cap activity triggered by primary electrons. The energy 
of electrons is only a few times higher than the thresh- 
old energy required to induce pair creation in the presence 
of a dipolar magnetic field, with other restrictions applied 
when necessary. Electromagnetic cascades induced via cur- 
vature radiation were treated in the same way as described 
by DH82. The important ingredient of the model is the as- 
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sumption that secondary particles, produced in cascades due 
to one-photon absorption, contribute to overall gamma-ray 
emission similarly as primary electrons. 

The model was confronted with the gamma-ray lumi- 
nosities for seven pulsars inferred from available data from 
OSSE , COMPTEL , and EGRET experiments. We find that 
the model is consistent with the existing data. Moreover, the 
model does not lead to any violence of energetics for pulsars 
with low spin-down luminosity L s( j , - the predicted gamma- 
ray luminosity L 7 never reaches L s( j. It avoids, therefore, the 
problem of 'an empirical gamma-ray death-line' as raised by 
Arons (1996). We also used the EGRET archive of 350 up- 
per limits along with OSSE , and COMPTEL upper limits 
(published for 15, and 18 pulsars, respectively), to find likely 
restrictions on the model. We have used some updates with 
respect to the data used by Nel et al.(1996) in their analysis. 
The EGRET upper limit for B1951+32 was replaced with 
its detections by EGRET (Ramanamurthy et al. 1995), and 
COMPTEL (Kuiperet al. 1996a). Similarly, the EGRET up- 
per limit for B 1509-58 was replaced with OSSE detection 
(Schroeder et al.1995). In the case of B1929+10 the model 
distance of 170 pc was replaced with 250 pc (see Yancopou- 
los, Hamilton & Helfand 1994 for detailed arguments). 

For a fixed value of L s d , the predicted L 7 depends rather 
weakly on magnetic field strength B as long as 10 11 G 
B <, 10 13 G (Dyks 1997), especially for pulsars with L sd ^ 
10 34 erg s _1 . That is why all evolutionary tracks calculated 
for high values of B converge roughly to the asymptotic re- 
lation of equation ^ for L s d ii 10 34 erg s _1 (dashed line in 
Fig. 3). Only as B enters the domain of millisecond pulsars 
(~ 10 8 — 10 9 G), L 1 drops significantly. Therefore, all mil- 
lisecond pulsars, including J0437-4715, are expected to be 
very weak gamma-ray emitters regardless of their L a( j, and 
their EGRET upper limits alone are still one order of magni- 
tude above our predictions for their gamma-ray luminosities. 
Qualitatively quite similar behaviour of millisecond pulsars, 
though for different physical reasons, results from the model 
of Dermer & Sturner 1994, which was used explicitely in the 
context of millisecond pulsars by Sturner & Dermer 1994. 
For their luminosity L smA = 1.1 x 10 10 £ 3/2 P~ 3 erg s" 1 of 
gamma-rays beamed into a solid angle f2 7 « 1.5 x 10 _3 P _1 sr 
(eqs.2 and 3 of Sturner & Dermer 1994, respectively), the 
apparent gamma-ray luminosity for 1 sr can be expressed as 

L 7 = fi- 1 L S £>94 « 1.2 x 10 9 B 1/2 L^ 2 erg S "\ (9) 

and accordingly for two pulsars with B = 10 9 G and 10 12 G 
but identical spin-down luminosity L s d the former object 
will be placed below the latter one in a diagram like Fig. 3. 
The EGRET upper limit for J0437-4715 (~ 10 32 erg s" 1 ) is 
well above 1.6 x 10 30 erg s _1 resulting from equation (^|). 

Out of several objects in the analysis of Nel et al.(1996), 
with uncomfortably low EGRET upper limits, and contra- 
dicting thus the models they discuss, only two are left as a 
potential threat to our model: B0950+08, and B0656+14. 
The EGRET limits for these two sources are too low to be 
accommodated by the model. It is encouraging, however, 
that B0656+14 was reported as a possible EGRET source 
(Ramanamurthy et al. 1996) . Both pulsars were on the prior- 
ity list of COMPTEL but with low ranks, and no results are 
available for them so far. There are no OSSE limits available 
for B0656+14 either. Definitely, B0656+14 deserves more 
attention as a promising target for gamma-ray experiments 



below the energy range of EGRET . Its parameters are very 
similar to those of Geminga and B1055-58. Moreover, all 
three pulsars are strong X-ray emitters, and are thought to 
be the best candidates for initial cooling (Ogelman 1995, 
Becker & Triimper 1997). On the other hand, the com- 
bined upper limits available for B1046-58, and B1929+10 
(EGRET , COMPTEL , OSSE in both cases) don't rule out 
our model. 

The upper limits for energy fluxes adopted from Nel 
et al.(1996), and used also in this work require a word of 
comment. They were inferred from upper limits for photon 
fluxes under an assumption that all photon spectra above 
100 MeV have a spectral index 7, obeying a trend derived 
by Thompson et al.(1994) from five EGRET pulsars: 

7 = 0.331ogr-3.08, (10) 

where the characteristic age of pulsars r = P/2P is ex- 
pressed in years. The trend is based essentially on the Crab 
pulsar (7 = 2.16, r = 1.3 x 10 3 ) on one side, and on B1055- 
52 (7 = 1.18, r = 5.3 x 10 5 ) on the other side. With the new 
determination of the spectral slope for B1055-58, 7 = 1.59 
(Fierro 1995), the prescription for 7(7") looks questionable. 
As a consequence, upper limits for L 1 derived for old pul- 
sars, especially for millisecond pulsars, might be somewhat 
tighter. That would put models discussed by Nel et al.(1996) 
into even deeper trouble, whereas the model we propose still 
remains intact. 

Pulsars from the database of Taylor, Manchester & 
Lyne (1993) extended by Taylor et al. (1995), plus Geminga, 
arranged in a traditional ranking, based just on spin-down 
fluxes L B <i/D 2 , start with Crab and five other gamma-ray 
pulsars. But then, there is a wide gap (of no gamma-ray 
detections) before the seventh gamma-ray pulsar, B1055-52, 
emerges as No. 33. The gap contains several millisecond pul- 
sars, with their flagship J0437-4715 taking very high overall 
position - No. 7. 

Our ranking of top 30 candidates for gamma-ray emis- 
sion, arranged by a predicted flux resulting from equations 
§ and (§), 

771 t 1/2 

h = c n± - E ^ L « , (11) 

is presented in Table 2. The pulsar database of Taylor et 
al.(1993, 1995) ordered by / 7 starts now with Vela, then 
goes Crab and Geminga. B 1055-52 advances by 13 positions 
to No. 20. The millisecond pulsars (from the gap), including 
J0437-4715, disappear from the list of 'Top 30'. 
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Table 2. Our ranking of top 30 gamma-ray candidates. The seven 
CGRO detections are marked with 7. 





PSR B 


PSR J 




D 
[kpc] 


L~f/D 2 

r 1 On 

[ erg s cm J 


1 


0833-45 


0835-4510 


/ 


0.50 


0.1302E-07 


2 


0531+21 


0534+2200 


'Y 


2.00 


0.8963E-08 


3 




0633+1746 


I 


0.15 


0.7308E-08 


4 


1706-44 


1709-4428 


I 


1.82 


0.6430E-09 


5 


1929+10 


1932+1059 




0.25 


0.5678E-09 


(5 


0950+08 


0953+0755 




0.12 


0.5646E-09 


7 


1951+32 


1952+3252 




2.50 


0.4108E-09 


8 


0656+14 


0659+1414 




0.76 


0.2130E-09 


9 


1509-58 


1513-5908 


T 


4.40 


0.2080E-09 


10 


1046-58 


1048-5832 




2.98 


0.1904E-09 


11 




2043+2740 




1.13 


0.1653E-09 


12 


1823-13 


1826-1334 




4.12 


0.1239E-09 


13 


1800-21 


1803-2137 




3.94 


0.1193E-09 


14 


0740-28 


0742-2822 




1.89 


0.1143E-09 


15 


0114+58 


0117+5914 




2.14 


0.1136E-09 


1(5 


1757-24 


1801-2451 




4.61 


0.9188E-10 


17 




1908+0734 




0.58 


0.9129E-10 


18 


1727-33 


1730-3350 




4.24 


0.7281E-10 


19 




0538+2817 




1.77 


0.6567E-10 


20 


1055-52 


1057-5226 


7 


1.53 


0.6346E-10 


21 


0823+26 


0826+2637 




0.38 


0.6209E-10 


22 




1918+1541 




0.68 


0.5512E-10 


23 


0355+54 


0358+5413 




2.07 


0.4791E-10 


24 


1853+01 


1856+0113 




3.30 


0.4525E-10 


25 


0450+55 


0454+5543 




0.79 


0.4274E-10 


2(5 


1822-09 


1825-0935 




1.01 


0.4079E-10 


27 


1702-19 


1705-1906 




1.18 


0.3806E-10 


28 


1133+16 


1136+1551 




0.27 


0.2123E-10 


29 


0906-17 


0908-1739 




0.63 


0.1433E-10 


30 


1451-68 


1456-6843 




0.45 


0.8666E-11 
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ADDENDUM 

Already after submitting our paper for publication Lucien 
Kuiper pointed to us that the COMPTEL group has found 
indications for a signal from B0656+14 in the energy interval 
of 10-30MeV (Hermsen, W., et al., 1997, Proceedings 2nd 
INTEGRAL Workshop, ESA SP-382, 287). 
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